April, 1998 
LPTENS-98/11 

LEVEL SPACING OF RANDOM MATRICES IN AN 

EXTERNAL SOURCE 

E. Brezin a ) and S. Hikami^ 



a ) Laboratoire de Physique Theorique, Ecole Normale Superieure 

24 rue Lhomond 75231, Paris Cedex 05, France [] 
^ Department of Pure and Applied Sciences, University of Tokyo 
Meguro-ku, Komaba, Tokyo 153, Japan 



Abstract 

In an earlier work we had considered a Gaussian ensemble of ran- 
dom matrices in the presence of a given external matrix source. The 
measure is no longer unitary invariant and the usual techniques based 
on orthogonal polynomials, or on the Coulomb gas representation, 
are not available. Nevertheless the n-point correlation functions are 
still given in terms of the determinant of a kernel, known through an 
explicit integral representation. This kernel is no longer symmetric 
though and is not readily accessible to standard methods. In partic- 
ular finding the level spacing probability is always a delicate problem 
in Fredholm theory, and we have to reconsider the problem within 
our model. We find a new class of universality for the level spacing 
distribution when the spectrum of the source is ajusted to produce 
a vanishing gap in the density of the state. The problem is solved 
through coupled non-linear differential equations, which turn out to 
form a Hamiltonian system. As a result we find that the level spacing 
probability p(s) behaves like exp[— Css] for large spacing s; this is 
consistent with the asymptotic behavior exp[— Cs 2/3+2 ], whenever the 
density of state behaves near the edge as p(X) ~ X 13 . 
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1 INTRODUCTION 



The level spacing distribution p(s), first discussed by Wigner |IJ for nuclear 
energy levels, has been studied extensively in random matrix theory || and 
the universality of p(s) has been discussed for many cases, including the the 
distribution of zeros in Riemann zeta function [Q. The calculation of this level 
spacing distribution is always much more delicate than that of the n-point 
correlation functions. Indeed when two neighboring levels are separated by 
some interval s, it implies that all the other eigenvalues are outside of this 
interval, and consequently it involves all the correlation functions of those 
eigenvalues. In the simplest Gaussian ensemble it took many years of devel- 
opment of the theory of Fredholm determinants, Dyson's inverse scattering 
approach, tau-functions , before this problem was understood. For the Airy 
kernel which appears for the spacing in the vicinity of the edge of Wigner's 
semi-circle, Tracy and Widom developed a technique through coupled non- 
linear differential equations, which we have generalized here and applied to 
our previous work on random matrices in a matrix source Q. 

The problem we had studied concerns a Hamiltonian H = H + V, in 
which Hq is an iV by iV non-random matrix with a known spectrum, but V 
is a random Gaussian potential. The probability measure for H, beeing then 
a Gaussian in H — Ho, is not unitary invariant. Consequently the standard 
approach, which consists of tracing out the unitary degrees of freedom, and 
then using orthogonal polynomials, or a Coulomb gas representation, is not 
available. However we have obtained exact formulae, i.e. valid for matrices of 
finite size, for the n-level correlation functions with the help of the Itzykson- 
Zuber formula ||. It is remarkable that these n-point functions are still given 
by a determinant of an n by n matrix whose elements are given by a kernel 
K(X, //), a well-known fact when the orthogonal polynomials are available. 
However in our problem, this kernel is no longer symmetric, but we know an 
exact integral representation for finite N. 

If we had considered the level spacing distribution for those ensembles 
with a source, around a regular point of the spectrum, the usual proofs 
of universality would apply and the end result, in the appropriate scaling 
limit, would be identical to the Wigner ensemble level spacing. However near 
singularities of the spectrum new universality classes appear. For instance 
near the edge of the Wigner semi-circle, in a region of size iV~ 2//3 , a new kernel 
expressed in terms of Airy functions controls the correlation functions; Tracy 
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and Widom have succeeded to find the level spacing for this new kernel. 

In the source problem, we have shown earlier that we can have gaps 
in the spectrum of H, when the randomness of V is not strong enough to 
bridge the gaps between the eigenstates of H . If we tune the randomness in 
order to reach the limiting situation in which a gap closes, a new universal 
singularity appears at which the density of eigenvalues p(A) vanishes like 
|A| 1/3 . In this paper we obtain the level spacing distribution by an application 
of Fredholm theory from which we derive non-linear differential equations, 
which remarkably form a Hamiltonian system. This generalizes earlier work, 
and makes it clear that the technique is general. For the Wigner case, with 
the sine-kernel obtained by Dyson, 

sin7r(a; — y) ,„ „ . 

K(x,y) = — r ± ^ , 1.1 

ir(x - y) 

Jimbo et al || had obtained long ago a closed equation for E(s), the prob- 
ability that the interval (— s/2, s/2) is empty. Tracy and Widom [0] had 
considered the Airy kernel 

K(x, y) = Mx)A l '(y)-A/(x)A l (y) 
x-y 

(where Ai(x) is an Airy function which satisfies Ai"(x) = xAi(x)), and ob- 
tained E(s) for the semi-infinite interval (s, oo) , when s is near a singular 
edge of the semi-circle. The method developed in this paper is easily appli- 
cable to those two cases, which are briefly reviewed in the Appendices A and 
B. 

In all those earlier cases, as well as for our source problem, the probability 
E(s) of emptiness of the interval (-s/2, s/2), is a Fredholm determinant, 



00 (— \) n rb rb 

E ( s ) = E / • ■ ■ / n£ =1 cfe fc det[i^,^)] M=1 ,..., n (1.3) 

n=0 nl Ja Ja 



if we choose (a = —s/2, b = +s/2) for the interval (a,b). The sine and Airy 
kernels are symmetric : K(x,y) = K(y,x). For the source problem K(x,y) 
is no longer symmetric under exchange of x and y, although its square under 
convolution yields back the kernel itself as for the simpler cases]!]]. Starting 
with Fredholm theory for E(s) and extending the analysis of Tracy and 



3 



Widom [0, §, ||, we derive for our problem a Hamiltonian system, leading 
to coupled non-linear differential equations. From those equations one can 
determine the function E(s), and in particular its asymptotic expansion for 
large s, which is interesting. For the sine kernel, it behaves as E(s) ~ 
exp(— \s 2 ); for the Airy kernel, it becomes E(s) ~ exp(— ^s 3 ). In our 
previous workQ, we had applied a Pade analysis to a function R(t) related 
to E(s), 

E(s) = expff R(s)ds] (1.4) 
Jo 

where the variable s is defined as 

rs/2 

s = p{x)dx (1.5) 

J-s/2 

We had made the ansatz in our previous paper [[| that R{s) behaves like s in 
the large s limit, namely that E(s) is Gaussian in terms of s. More generally 
the ansatz was that for a density of state behaving as p(x) ~ x@, then 

E(s) ~exp[-Cs 2/3+2 ] (1.6) 

Indeed this ansatz agrees with the results for the sine-kernel, for which f3 = 0, 
and for Airy kernel (3 = 1/2. It also agrees with our gap closing singularity 
for which the density of state gives an exponent (3 = 1/3. The result of the 
present analysis confirms this conjecture since it gives E(s) ~ exp[— Cs 8 ^ 3 ]. 



2 The kernel at the closure of a gap 

We consider an iV x Af Hamiltonian matrix H = Hq + V, where Hq is a 
given non-random Hermitian matrix, and V is a random Gaussian Hermitian 
matrix 0, [H], 12, 13], 14]. The probability distribution P(H) is thus given by 

P(H) = le-*™* 

= _L e -fW-2«,H) (21) 

z 

For the n-point correlation functions, defined as 

i?n(Ai, A 2 , ■ • • , A n ) =< ^Tr5(Ax - H)^Ti5(X 2 - H) ■ ■ ■ ^Tr5(A n — H) > 

(2.2) 
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we have derived earlier [|1^, |TJJ] the expression 



i2 2 (A l5 A 2 ) = K N {X lt \i)K N (\ 2 , A 2 ) - K N (X lt X 2 )K N {X 2 , X 1 ) (2.3) 
with the kernel 



^(Ai,A a ) = / - I — TT(^— -)— e-f- 2 -^-^^ 

J 2,71 J 2ni zs x u — a 7 u — it 

(2.4) 

Similarly the n-point functions are given in terms of the determinant of the 



nxn matrix whose elements are given by this same kernel Kw(Xi, Xj) [15 
In this kernel Kn(Xi, A 2 ) was considered in the large N limit, for fixed 
N(Xi — A 2 ). In this limit one can evaluate the integrals Q2.4j) by the saddle- 
point method. The result was found to be, up to a phase factor that we omit 
here, 

^tv(Ai, A 2 ) = -— sin[7ryp(A 1 )] (2.5) 

where y = N(Xi — A 2 ). Apart from the scale dependence provided by the 
density of state p, the n-point correlation functions have thus a universal 
scaling limit, i.e. independent of the deterministic part H of the random 
Hamiltonian. 

In order to generate a tunable gap we consider the simple case for which 
the eigenvalues of Hq are ±a, each value being N/2 times degenerate. When 
the randomness of V is small, the average density of eigenvalues is made of 
two disjoint segments located around the points ±a. When the randomness 
increases one reaches a critical point at which the gap closes. When that 
happens in the vicinity of the origin of size iV~ 3//4 we find a new class of 
universality for the density of states and for the n-point functions^]. Then 
the kernel (12.41) becomes 



K N (X lt X 2 ) = (-1)- +1 / /^(4 ± 4) f J—e-^-^-^^ 

J In J ztti u z — a z u — it 

(2.6) 

From this expression, we obtain the density of state p(X) = K N (X,X). The 
derivative of p(X) with respect to A eliminates the factor u — it in the de- 
nominator of (|2.4|), and leads to the factorized expression, 



5 



where 



jf-jftpw = -mm (2.7) 

0(A) = r ^e-^fHa^-mtx (2g) 



Vi(A) = / — e-f « 2 ~f m(- a -« 2 )+^ (2.9) 

J 27TZ 

For N large, we may apply to the two integrals defining the functions <fi and 
ip the saddle-point method. When Ai and A2 are near the origin the saddle- 
points in the variables t and u move to the origin. Therefore for obtaining 
the large N behavior of near A = we can expand the logarithmic term 
in powers of t. One sees readily that the coefficient of the quadratic term in 
t 2 vanishes for a = 1; in fact three saddle-points are merging at the origin 
when a reaches one. This is the critical value at which the gap closes. We 
must then expand in the exponential up to order t 4 and we obtain 



^ dt 

1 

to 27T 



0(A) = f ^L e -^-mtx (21Q) 



Rescaling t to N l l 4 t' and setting A = N 3//4 x we find that 

0( x ) = N^^N-^x) (2.11) 

has a large N, finite x, limit given by 

i r°° 1 4 

0(x) = - / dte~^ cos(tx) (2.12) 

7T JO 

It is immediate to verify that it satisfies the differential equation, 

j/"(x) = x(f)(x) (2.13) 



From the integral representation (|2.8|) we obtain easily the Taylor expansion 
of this function at the origin 



y/2 « r(| + f)(-i) m 2 
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and its asymptotic behavior at large x, 

<l>(x)~d—x 3e cos (__ a;3 - ) (2.15) 

For the second function (2.9|), in the scaling limit, N large, A small, iV 3 / 4 A 
finite, we may expand up to order u 4 , and define 

^(x) = N 1/4 ip(N- 3/4 x). (2.16) 

In the large N, finite x, limit we find 

$( x ) = I ^eT+». (2.17) 

Jc 27TZ 

The integral is over a path consisting of four lines of steepest descent in 
the complex u-plane. Along these straight lines, the variable u is changed 
successively into e <'u and e ~ l u. This leads to 



A /^J />00 4 

ij}(x) = -Im[- / e -^(e x ^ - e- XUUJ )} (2.18) 

7T JO 

in which lo = e~i . The function ty{x) satisfies the differential equation, 

fy"{x) = -xip(x) (2.19) 



and again we find from (|3.19|) the Taylor expansion 



1 00 (-D n r 4n+1 (2nV 
^( x ) = _^vL_^ tP 1 (2-20) 



and the large x behavior 

^ 2l /A z -| e M cos(^xt + — ) (2.21) 
V 37r 8 3 

As shown in ||, one may express the whole kernel K N (Xi, A 2 ) of ( |3.1| ) in 
terms of the two functions <j) and ■?/> in the scaling limit. Defining 

\ 1 = N-V 4 x, A 2 = N-*' A y (2.22) 
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K(x, y) = N^K^N^X!, iV 3 / 4 A 2 ). (2.23) 
in the large N, finite x and y, limit, we had shown that 

x-y 

Note that if follows from ( |2.12j ) that <p{x) is an even function, whereas ip(x) 
is odd ( |2.18D . It implies in particular 



K(-x,-y)=K(x,y). (2.25) 

Therefore the density of state is given by 

p( x ) = -[^( x )ip"(x) - 0"(x)V/(x) + x4>(x)i>{x)} (2.26) 

In the large x limit, it behaves as p(x) ~ x 1//3 . Hereafter, we denote simply 
and if) by and ip, respectively. 



3 Fredholm theory 

The level spacing function E(s), the probability that there is no eigenvalue 
inside the interval (-s/2, s/2) centered around the singular point s = 0, is 
given by the Fredholm determinant, 

E(a,b) = det[l-K] 

00 ( — \) n rb rb 

= ■■■/ n« =1 rfx fc det[K(x 4 ,x,)] M=1 ,..., n (3.1) 

n=0 n[ Ja Ja 

if we choose for (a,b) the interval (-s/2, s/2). The sine and Airy kernels 
are symmetric kernels, K(x, y) = K(y, x), whereas our kernel in ( [2.24|) is not 
symmetric, since the two functions 4>(x) and ip(x) are different. Our notations 
here are as follows : the kernel K(x,y) is defined by (p. 24 ). However it acts 



on the interval (a,b) ; therefore we have used in ( |3.1| ) the kernel K(x, y) 
defined by the restriction of K to the interval : 

K(x, y) = K(x, y)9(y - a)6(b - y) = K{x, y)Q(y) (3.2) 
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in which 6{x) is the Heaviside function, and we have used for convenience 
the notation 

Q{y)=9{y-a)9{b-y). (3.3) 

In order to calculate the derivative of the logarithm of E(a, b) with respect 
to the end points one writes 

\nE(a,b) = Trln(l — K), (3.4) 

and thus 

db ~ 1-K0b" 1 ' 

^From (|3.2|) we find 

?^l = K(x,b)S{y-b), (3.6) 

and therefore, if we introduce the Fredholm resolvent K (6, b) , defined by 

K , , 

K = ^ 3.7 

l-K 

we obtain 

^f^ = -A-(M) (3.8) 
Similarly the derivative with respect to a is 

*!§M = *(.,.) (3.9) 

Therefore, when we choose for (a,b) the interval (-s/2,s/2), we have 
dhiE(s) l.d <9 , , s, ~,s s. 

This leads to 

E(s)= exp[- f°K^, S -)ds>] (3.11) 
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We now define six functions, obtained by acting with the operator (1 — 
K)^ 1 on the functions <f>(x), ip(x) and their first two derivatives, 

q (x) = (1 - K)- 1 ^) (3.12) 
= <j>(x)+ f K(x,y)<f)(y)dy+ [ f K(x,y)K(y, z)<f)(z)dydz + ■ ■ ■ 

J a J a J a 

We find it easier to use Dirac's notations 



q (b,a;x) =< x\ 1 f \<j> > . (3.13) 
1 — K 



Similarly we define 



q n (b,a;x) =< x\ — > (3.14) 

where (fi( n \x) is the n-th derivative of (f>(x). For the function ip(x), we act 
with the operator (1 — K)~ l on bras, i.e. dual vectors, rather than on kets. 
Because of the lack of symmetry of the kernel K this introduces a kernel L, 
and 

p n (b,a;x) = (-1)™- 1 < ^ 2 - n )\-^\ x > (3.15) 

1 — L 

where 

L(y,x) = Q(y)K(y,x) (3.16) 

When we set x = b, and a = —b, then the six functions q n and p n become 
functions of the single variable b, which we denote as 

Qn(b) = q n (b, -b- 6), P n (b) = Pn (b, -b; b) (3.17) 

The calculation of the derivatives of these functions Q n {b) and P n (b) im- 
plies to consider separate variations in the functions q n (b, a; x) and p n {b, a; x) 
with respect to b, a, and x, before we set a = —b and x = b. The calculation 
is tedious, but not difficult, and we have given more details in the Appendix 
C. The resulting differential equations are 

2 

Qo — Qi + tQiPiQo 
b 

Qi = Q2-^QiPi-Qou 
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Q 2 
Pi 

P-2 



bQo + TQ1P1Q2 - Qiv 
b 

-bP 2 - IQ1P1P0 + Pm 
b 



(3.18) 



QiP? + P2V 



-Pi - tQiPiP* 



where a dot means taking the derivative with respect to b. The two auxiliary 
functions u and v are defined as 



u =< tp\qi >, v =< ip'\qo > 

and they satisfy 

u = -2P 2 Q U v = 2P 1 Q . 
From these equations, we find succesively the relations 

u + v = -2P 2 Q , 

then 



Qo = Qi(i + ^ 



p 2 = - ¥ 



Using ( |3.20|) and (|3.21| ), we obtain from the second equation of ( |3.18|) , 

2 vu 2 



2P2Q2 



an 



u 



u + v bu + v 
and using the fifth equation of Q3.18D , we have 

2 iiv 2 



u{u + v) 



2Q P = v 



uv 



u + v bu + v 



+ v(u + v). 



(3.19) 

(3.20) 

(3.21) 
(3.22) 

(3.23) 
(3.24) 



Taking the derivatives of these two equations and using the third and sixth 
equations of fl3.18|) , we end up with two coupled non-linear equations 



d 3 u ,2u 1 
— + (— -l)[b( u + v) + —— 1 
db A b u + v 



vu + 2vu) 



uv fn . ... 2v(u) 2 
(2v + u)}-' 



(u + v) 



b 2 (u + v) 







(3.25) 
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d 3 v ,2v , ri . , 1 , uv , 2w(t>) 2 

— + — + 1 ) 6 (m+k — (uv + 2uv) + - (2u + v + ^ =0 

db A b u + v (u + v) 2 b 2 (u + v) 

(3.26) 

In the large 6 limit, the asymptotic behavior of the solutions of these equa- 
tions is obtained under the form 

U =J + ^ 2/3 + " " " (3-27) 
b 2 A 

v = -- + -b 2/3 + ■■■ (3.28) 

4 2 v ' 

Inserting these expressions in ( |3.25| ), we find A from the coefficient of the 

terms of order b 1 ^ 3 , 

A = ~(\) 1/3 (3-29) 

The kernel K(b, b) may then expressed as 

K(b,b) = bP 2 Q + Q 2 P 1 + Q 1 P -uP 1 Q -vP 2 Q 1 

- ^(ifQi-^-QoA) 2 (3.30) 

Remarkably this kernel at coinciding points is in fact a Hamiltonian, that we 
denote as 

H(b) = K(b,b), (3.31) 
from which Hamilton's equations 

^ = W n ' = ~Wn (3 ' 32) 

coincide with the differential equations ( |3.18| ). This allows one to obtain a 
simple expression for its derivative with respect to b becomes 

dH(b) dH 
db db 



Q P 2 + ^P 2 Qj (3.33) 



In terms of u and v , it becomes 



dHjb) u + v (uv) 2 
db 2 2b 2 (u + v) 2 { ' 
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Thus, from the previous result, we get the large b behavior 

dH(b) 



db 



5 • 2- n / 3 6 2 / 3 (3.35) 



. Integrating once we find for the hamiltonian H(b) ~ 3 • 2~ n / 3 5 5 / 3 ~ 3 ■ 
2-16/3 s 5/3 j n i ar g e s ii m it. 

From (|3.11|) , we have 

E(a) = expf- f H(s')ds'] 
Jo 

~ exp[-9-2- 25 / 3 s 8 / 3 ] (3.36) 

We have thus derived the exponent s 8//3 which was a mere conjecture in 
our previous article jp. There we had performed a simple Pade analysis of 
the small s expansion, and assumed that it was Gaussian at large s in the 
variable s. This had led us to the estimate 

E(s) ~ exp[-0.332S 2 ], (3.37) 

or since s ~ ^(^) 4/3 s 4/3 , E(s) ~ exp[-0.0358s 8/3 ]. Our analytic result 
( |3.22| ) gives E(s) ~ expf— 0.0280s 8 / 3 ]. Thus the estimation by a simple Pade 
analysis was not too far from the exact result. 



4 SUMMARY AND DISCUSSION 

In this paper, we have investigated the level spacing probability for the case, 
in which two edge singularities collapse. By the use of Fredholm theory, we 
have derived an expression for the level spacing probability, whose logarith- 
mic derivative turns out to act as a Hamiltonian. The same strategy allows 
one to solve as well the simpler cases, such as the sine-kernel (relevant to 
the level spacing for ordinary non-singular) points of the spectrum and the 
Airy kernel (which applies to a single edge singularity). The corresponding 
Hamilton's equations determine fully the level spacing, and in particular one 
can obtain analytically its asymptotic expansion at large spacing s. This 
allowed us to confirm the conjecture that we had made in our previous arti- 
cle, on the asymptotic Gaussian behavior of the level spacing in terms of the 
variable s = fl% 2 p(x)dx, which is the number of eigenvalues in the interval 
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of size s J3J. We have thus derived here that the level spacing probability 
E(s) behaves exp[— Css], with a constant C that we have analytically de- 
termined. More generally the three cases that we have solved are consistent 
with the asymptotic Gaussian behavior of E(s) with respect to s, i.e. to 
E(s) ~ exp[— Cs 2/3+2 } for large spacing s, whenever the interval of size s is 
around a point at which the density of state behaves as p(A) ~ X^. The 
scaling variable s is in fact related to the actual interval between the eigen- 
values by absorbing a power N^^ +1 \ If we let N go to infinity first, at 
fixed interval, in which case s is also large, a finite limit of ln ^^ for large 
N, implies that E(s) does fall for large s as E(s) ~ exp[— Cs 2l3+2 ]. In other 
words, E(s) behaves for largeiV, fixed interval, as a partition function. 
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Appendix A: Level spacing probability for the sine kernel 

The sine kernel is defined by 

x-y 

where 4>(x) = sinx satifies <p"{x) = —<p(x). (For convenience we have ab- 
sorbed in the normalization the usual factor ir). We consider q(x) and p(x) 
defined by 

q(x)=<x\ T> \4>>, p(x) =< <f>'\ - 1 f \x > (A.2) 
1 — K 1 — K 



and the Fredholm resolvent K 



We have 



K 

K = (A.3) 

l-K 



(x-y)K = <x\[X,k\\y>=<x\[X,—Lr]\y> 

1 — K 

= <x\ — —r\X,It\^—r\y> (A.4) 

The definition (|A.1|) of the kernel reads 

[X,K\ = \(f>><(l/\-\<f/><<f>\ (A.5) 

since 

(x-y)K(x,y)=<x\[X,K]\y> (A.6) 
Thus, from (|A.4j ), we obtain 

K( x ,y) = q W p M-«<»M x ) (A.7) 
x-y 

Since the functions q(x) and p(x) depend upon the interval (a,b), we 
denote them more precisely as q(b, a; x) and p(b, a; x). We then set x = b and 
vary b, i.e. take the derivative of q(b, a; b) at fixed a 

dq(b, a:b) , . „ 1 , , 1 dK 1 . , . . „. 

L =< b \ D - 0> + <6 A - 8 

db l-K l-K db l-K 
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where D is the derivative operator: < x\D\f >= f'(x). From the definition 
of K we have, 



dK(x,y) 
db 



K{x,y)5{b-y) = K\b >< b\ 



(A.9) 



Thus the second term of ( |A.8 ) becomes K(b,b)q(b). The first term of ( |A.8| ) 
becomes 



1 



<b\D -|0 > = p(b)+<b\[D, 



1-K 



1-K 



> 



p (b)+ < b\^-Lr[D,k]^ 1 r\ ( f>> (A.10) 



Since 



<x\[D,K]\y> = (— + —)K(x,y) 



K(x,y)(5(y-a)-5(y-b)) 



we obtain 



[D, K] = K\a>< a\- K\b >< b\ 
/From ([Q| ) and flATTOP , we have 

dq(b, a; b) 



db 

Similarly, we have for p(b) as 

dp(b, a; b) 
db = 



p(b, a; b) + K(b, a)q(b, a; a) 



-q(b, a; b) + K(b, a)p(b, a; a) 



We denote q(b, —b; b) by Q(b), 

dQ(b) dq(b, a; b) 



db 



db 



a=—b 



dq(b, a; b) 
da 



a=—b 



Since 



dq(b, a; b) 
da 



< b -(-5-) ^ > 

—K(b, a)q(b, a; a) 



(A.ll) 
(A.12) 

(A.13) 
(A. 14) 
(A.15) 



(A.16) 
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we have 

Q(b)=P(b) + 2K(b,-b)Q(-b) (A.17) 

where 

K (b ,- b)= 9mw (A . 18) 

b 

Note that Q(—b) = —Q{b) and P(—b) = P(b). Similarly we have an equation 
for P(b). Thus we obtain finally 

q = m- 2 -^) 

2P 2 

P = Qi-^-l) (A.19) 
The function K(b, b) is related to P and Q as 

K(b,b) = P 2 + Q 2 -^-^- (A.20) 

which gives the logartihmic derivative of the level spacing probability E(s). 
Noting b = s/2, we have 

f ■ f^-» 

In the large s limit, we have Q ~ P ~ s 1 / 2 /2. From ( A.20Q , we obtain 
H(s) = K(b,b) ~ s/4. The level spacing probability E(s) behaves thus in 
the large s limit as 

E(s) ~ exp[- / -ds'l 



o 4 

exp[-^s 2 ] (A.22) 



For small s, we have by solving ( |A.21| ) iteratively, with the initial condition 
P(0) = 1, Q(0) = , 

« = | - ^ + 0(s 5 ), P=l + s+ 7 -s 2 + 0(s 3 ) (A.23) 
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Then from ( |A.20| ), H (b) = K(b, b) behaves for small s as 

H(s) = 1 + s + s 2 + 0(s 3 ) (A.24) 

This leads to 

E(s) = expf- f H(s')ds'} 
Jo 

= l-s + 0(s 4 ) (A.25) 

which is consistent with all the well-known results on this well-studied case 
[0]- 

Appendix B: Level spacing probability for the Airy kernel 

The sine kernel applies to a regular point of the spectrum. However 
when one studies the vicinity of the edge of the spectrum, in the appropriate 
scaling limit, the correlation functions and the level spacing are given by 
an Airy kernel. The level spacing has been studied by Tracy and Widom 
[0]. We repeat here the same technique. Consider the interval (— s/2,s/2) 
as in the sine case. We denote the Airy function Ai(x) by <p(x), it satisfies 
(j)"{x) = x<fi(x). We use the same notation q(x) and p(x) as in ( |A.2| ). As for 
the sine case, the Fredholm resolvent K (a, b) is given by 

K{aM = q[a)m ~ P } a)q{h) (B.l) 



We have 



(-^ + -^)K(x,y) = -<P(x)<P(y) (B.2) 



Therefore, we have 

[D, K] = — 10 >< 0|6 + K\a >< a\ -K\b>< b\ (B.3) 
Thus, by the same procedure as sine case, we get 

^=p(b)-qu + K(b,a)q(a) (B.4) 

^ = bq(b) + up(b) - 2q(b)v + K(b, a)p(a) (B.5) 
ob 
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where 

u =< <p\q >, v =< <p\p > (B.6) 
The Fredholm resolvent K(b, b) becomes 

K(b, b) = p 2 (b) - bq 2 (b) - 2up(b)q{b) + 2q 2 {b)v 
+^[q(b)p(a) - p(b)q(a)][q(a)p(b) - p(a)q(b)] (B.7) 

Again H(b) = K(b, b) acts as a Hamiltonian, since we have 

dH(b) _ Jqjb) dHjb) dp(b) 
dp(b) db ' dq(b) " 8b 1 ' 

The derivative of the Hamiltonian becomes 

dH(b) 2 (q(b)p(a) - p{b)q{a)f 

-dT = - q {b) + {b^ (B ' 9) 

In the Airy kernel, due to the parity around the edge, H(a) and H(b) are 
different. The quantity q(b) and p(b) becomes exponentially small in the 
large b limit (b > 0) as same as Airy function. We have for a = —s/2, 

9q ^ -- p{a) - q{a)u - K{a, b)q{b) (B.10) 



da 

dp(a) 



aq(a) + upia) — 2q(a)v — K(a, b)p(b) (B-H) 
da 

The Hamiltonian H(a) becomes 

H(a) = p 2 (a) — aq 2 (a) — 2up(a)q(a) + 2q 2 (a)v 

+ K(a,b)[q(a)p(b)-p(a)q(b)} (B.12) 

Since K(a, b) can be neglected for the large b limit, and we have a relation, 
u 2 — 2v = q 2 , we get 

= aq(a) + 2q%a) (B.13) 

aa A 

This leads to q (a) ~ y/s/2, and 

H(a) ~ s 2 /16 in the large s limit. Then, we get E(s) ~ exp[— s 3 /96]. In 
order to check the usefulness of the Pade analysis we expand now at small s. 
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The Airy function A(x) has the Taylor expansion for small x, 

1 4 

A(x) = c 1 [l + -a; 3 + -x 6 + ---] 
6 6! 

- c 2 [x + ^x' + ^x 7 + •■■] (B.14) 

where c x = 3~ 2 / 3 /r(2/3) and c 2 = 3 _1 / 3 /r(l/3). If we use 5 , defined in 
( |1.5| ), it is related to s in this Airy case by 

~ s = els + £l^ s 3 _ ^L s 5 _ _J^ S 7 + G( 9) (B 15) 

2 12 960 16128 v ; v ; 

The level spacing -E(s) is thus expanded in powers of s as 

E(s) = + + 

V ; 6 V 3 8 ' c\ 

= l-s + 0.544868s 4 -42.5418s 6 + O(s 8 ) (B.16) 
We now apply a Pade analysis to H(s), or rather to 

R(s) = ^-\nE(S) 
as 



1 + ai§ + a 2 s + a 3 s 
1 + bis + b 2 s 



(B.17) 



where we have a\ = 754.156, a 2 = 1640.76, a 3 = 1638.58, b\ = 753.156, and 
b 2 = 886.601. From E(s), differentiating twice with respect to s , we obtain 
p(s), which is slightly different from the usual 'Wigner surmise' function of 
the sine case. 

Appendix C: Level spacing probability for the gap closure ker- 
nel 

We consider now the kernel 

K , y) = (ai) 

x-y 

or, in operator notations, 

[X, K] = >< - |0" >< VI - |0 >< ip"\ (C.2) 



20 



Similarly to ( |A.4j ), we have 



(x-y)K = < x\—^[X,k]—^\y > 

= Qi(x)pi(y) - q 2 (x)p (y) - qo(x)p 2 (y) (C.3) 
The derivative of q n (b, a; b) for a fixed a becomes 

d -^f^ = < W i |^.) > + < t| * (f) i > 

96 1 - K \-K oh \-K 

= qn+l (b,a;b)+ Kbl—^ID^}—^^ > 

1 — K 1 — K 

+ <b\^^\b><b\—^\<p {n) > (C.4) 
l-K l-K 

We have also 

d d 

<x\[D,K]\y> = (— + —)K(x,y)+<x\K\aXa\y> 

- < x\K\b >< b\y > (C.5) 

The first term is simply —<p(x)ip(y). This leads to 

[D, K] = -|0 >< ip\Q + K\a >< a\ - K\b >< b\ (C.6) 

Therefore 

dg n (b,a;b) 
db 

and 



q n+1 + K(b, a)q n (a) - q (b) < i/j\q n > (C.7) 



q 3 = <x\^-^X\<j)>=<x\X^—^\(j)> + <x\[^—^,X]\<f>> 

= xq (x)+ < x\ — —,r\k,X] — > 
l-K l-K 

= xq - v 2 qi + Uiq 2 + v 3 q (C.8) 
where u\ =< i/;\q >, v 2 =< if)'\qo > and v 3 =< ip"\qo >■ 
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The function p n (x) is denned by 

p n (x) = (-IT' 1 < ^ (2 - n) \—^r\x > (C.9) 

1 — L 

where 

L(y,x)=G(y)K(y,x) (CIO) 

We have 

[D, L] = -Q\(f) >< + \a >< a\K -\bxb\K, (C.ll) 
in which 6 is a local operator defined by 

< y\G\y' >= 5(y - y')9(y - a)6(b - y) (C.12) 

Thus we obtain 

"' ] = -Pn-M-p^b) < V> (2 - n) ko > + Pn (a)K(a,b) (C.13) 



db 

with p-i(x) obtained as 

P-i(x) = -<ip"'\ — l —\x > 

1 — L 

i - i 

= xp 2 (x)— < — j[x,L]- — j\x> (C14) 
= -xpo(x) - pi(x) < V'lgi > -P2 < Vl<?2 > -Po(x) <ip\q > 

where 

<y\[X,L]\x> = (y-x)Q(y)K(y,x) 

= 9(y)(|0' >< v'l - 10" >< VI - 10 >< V"I)(C15) 

The function (f>(x) is an even function of x, and qo(x) becomes an even 
function. The function ip(x) is a odd function. Therefore, we have u\ = 
v 3 = for the interval (-b,b). Also we < ip\q2 >=< ip\qo >=< ip"\qo >= 0. 
Nonvanishing quantities are w 2 =< V|?i >> an d ^2 =< V'ko >• We denote 
them simply by u and v. 

Noting that 

n dqn \ d( * n \ (ria\ 
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we obtain ( |3.18| ). The equations for P n in ( |3.18| ) are also obtained similarly. 
The derivative of u = u 2 becomes 

u=-2P % {b)Q x (b) (C.17) 

This is obtained as 

t =< * |e T^§r^' >+ * (%i(6) (ai8) 

Using 

^ = K(x, y)5(y - b) = K\b >< b\ (C.19) 

we have 

- = - P2 (b) gi (b) (C.20) 
The derivative of u by a becomes, similarly 

^ = P2(o)gi(o) (C.21) 



da 

Putting a = —b, we have 



du du 

^ r%7 \a=—b ~7{ \a=—b 

ob oa 
2P 2 (b)Q 1 (b) (C.22) 



This is Eq.flngp. 

Appendix D: Modified kernel 

We have considered the case of the fixed external source eigenvalues at 
a = ±1. Here, we take this external eigenvalue a as a = 1 + 2N~ 1 ' 2 a. The 
parameter a measures the approach to the limit a = ±1 in the large N limit, 
where the gap is closed. This change of a modifies the function and 
i)(x). We have in the large scaling limit 

^(A) = r ^ e -£-crf»+«A (d.i) 
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and it satisfies 

0'" - 2a0' - X(f) = (D.2) 
The equation, which t/>(A) satisfies, is also modified as 

- 2aip' + Xtp = (D.3) 

Following the same procedure in the Appendix B in we have a kernel, 
which is slightly different from (|3.14 ), 



a;-y 

However, the density of state p{x) is expressed by the same equation as 
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p(x) = -[</>' (xW{x) - (f>"(x)i)'{x) + x(j)(x)ip(x)\ (D.5) 

The large A behavior of 0(A), for a fixed a, is obtained by a saddle point 
method. If we make a change of t by A 1 / 3 /!;, we find that the new term at 2 
becomes negligible compared with other terms, which becomes order of A 4 / 3 . 
Then, we obtain the large x behavior of p(x) as x 1 ^ 3 same as before. 

Using the same definitions for q n and p n , we have an equation, 

£>/ ,x + Qo(a)Po(b) + g 2 (a)P2(&) - 2aq (a)p 2 (b) 

K(a,b) = (D.6) 

a — o 

Although there are modifications in the differential equations for q(b) and 
p(b), the derivative of the Hamiltonian is given by ( |3.33|) , which gives the 



same asymptotic behavior of E(s). We note that the function <f>(x) in ( |D.1| ) 
appears for the second Painleve A 2 Gamier system (Appendix A in [[|), and 
the function <f)(x) satisfies the coupled linear partial differential equations 
about x and a. 
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